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Within the framework of a generalized factorization, semi-inclusive deeply inelastic scattering 
(SIDIS) cross sections can be expressed as a series of products of collinear hard parts and transverse- 
momentum-dependent (TMD) parton distributions and correlations. The azimuthal asymmetry 
(cos <f>) of unpolarized SIDIS in the small transverse momentum region will depend on both twist-2 
C^l , and 3 TMD quark distributions in target nucleons or nuclei. Nuclear broadening of these twist-2 

< — i . and 3 quark distributions due to final-state multiple scattering in nuclei is investigated and the 

' nuclear dependence of the azimuthal asymmetry (cos (f>) is studied. It is shown that the azimuthal 

asymmetry is suppressed by multiple parton scattering and the transverse momentum dependence of 
the suppression depends on the relative shape of the twist-2 and 3 quark distributions in the nucleon. 
A Gaussian ansatz for TMD twist-2 and 3 quark distributions in nucleon is used to demonstrate 
the nuclear dependence of the azimuthal asymmetry and to estimate the smearing effect due to 
fragmentation. 

Q_l' PACS numbers: 25.75.-q, 13.88.+C, 12.38.Mh, 25.75.Nq 
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I. INTRODUCTION 



In high-energy nuclear collisions, from deeply inelastic scattering (DIS) off a nuclear target to hadron-nucleus and 
■ nucleus-nucleus collisions, multiple parton scattering plays an important role in the reaction dynamics and the final 
' hadron spectra. It not only causes the transverse momentum broadening of the propagating parton but also leads to 
parton energy loss due to induced gluon brcmsstrahlung. A direct consequence of the transverse momentum broadening 
is the so-called Cronin effect or the broadening of the final hadron spectra in transverse momentum in both DIS off 
nuclear targets [l|, and hadron-nucleus collisions^ Q . Gluon bremsstrahlung induced by multiple parton scattering 
on the other hand leads to parton energy loss |5H8| and medium modification of the parton fragmentation functions 
a phenomenon known as jet quenching which is observed as the suppression of the leading hadron yields from 
parton fragmentation. The radiative parton energy loss and medium modification to parton fragmentation functions 
are also determined by the jet transport parameter or the average squared transverse momentum broadening per unit 
distance experienced by the propagating parton, which is related to the gluonic structure of the scattering centers 
in the medium 0, [uj. Therefore, extraction of the jet transport parameter through experimental measurements 
of transverse momentum broadening and jet quenching can provide important information about the properties of 
the medium, either cold nuclei or hot and dense QCD matter, as probed by energetic partons. The jet quenching 
phenomenon has been seen in high- energy heavy-ion collisions via strong suppression of not only the large transverse 
momentum single hadron spectra [H, but also the back-to-back dihadron correlation [lj]. The extracted jet 
transport parameter is found to be much larger than that in large cold nuclei as obtained by phenomenological studies 
of jet quenching in DIS flH [lfjj , indicating much higher gluonic density in the initial stage of high-energy heavy- ion 
collisions. 

Much efforts have also been devoted to the study of transverse momentum broadening in DIS off nuclear targets 
H, H, [l(J EE EE and the Drell-Yan dilepton production in proton-nucleus collisions [25rl27| within different 

approaches such as the color dipole modeljl7 Hl9fl and higher-twist expansion in the generalized collinear factorization 
formalism [2(J[2l[. It has been shown recently |24{ that the gauge invariant transverse momentum dependent (TMD) 
quark distributions in nucleons and nuclei can be expressed as a sum of higher- twist collinear parton matrix elements, 

tf(x,k±) = J ^e^ + y-(A | ^(O)^^-)-^^-) | A)5^(k±), (1) 
in terms of the parton transport operator 

W ± (y~) = iD ± (y-)+g f dCF +± (C), (2) 



where D± (y ) = d± + igA± (y ) is the covariant derivative. For brevity of presentation, we have used the light-cone 
gauge and all the transverse coordinates have been set to y±_ = 0± in the above. Since the above expression is valid 
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for both nucleon and nuclear targets, the nuclear dependence of TMD parton distribution functions will come from 
the nuclear dependence of the higher-twist parton matrix elements. Under the "maximal two gluon approximation" , 
the TMD quark distribution in a nucleus can be expressed as a convolution of the TMD quark distribution fq{x, £±) 
in nucleon and a Gaussian broadening [24| , 



The broadening width A 2 _f or the total average squared transverse momentum broadening, 

A 2F = / dZ N q F {Z N ). (4) 



is given by the quark transport parameter, 

<M6v) = ^/4(6v)[x/f WWo, (5) 

where Pn({;n) is the spatial nucleon number density inside the nucleus and fg (x) is the gluon distribution function 
in a nucleon, 

xfg{x) =- J 7^ eiXP+r ( N I ^(O)V(D I N). (6) 

A summation over the gluon's color index is implied in the above definition of the gluon distribution function. 

One can generalize the above approach to nuclear modification of higher-twist TMD parton distributions. In this 
paper, we will study the case of twist-3 TMD quark distributions which determines the azimuthal asymmetry (cos </>) 
of the unpolarized semi-inclusive deep-inelastic scattering (SIDIS) cross section defined with respect to the leptonic 
plane. Such asymmetry at large transverse momentum arises predominately from hard gluon bremsstrahlung in 
perturbative QCD [28|, [29| . Similar processes are also responsible for the azimuthal angle dependence of the Drell-Yan 
dilcpton production cross sections |30| and multiple parton scattering in p + A collisions can lead to a twist-4 nuclear 
dependence [3lT l32l]. In the small transverse momentum k± < 1 GcV/c region, however, the asymmetry is shown to 
arise mainly from the intrinsic transverse momentum of the quarks in the target nucleon or nucleus [33l [34| and is 
related to the higher-twist TMD parton distribution functions (35l - [38j . We will focus on the twist-3 contribution to 
the azimuthal dependence, proportional to (k±/Q) cos cj), in this paper. Here Q 2 = —q 2 and q is the four momentum 
transfer from lepton to target. Since the coefficient of the azimuthal asymmetry is related to twist-3 TMD parton 
distributions, one can take the same approach as in the case of the twist-2 TMD parton distribution to study the 
nuclear dependence of the twist-3 contribution to the azimuthal asymmetry of SIDIS cross sections. 

The remainder of this paper is organized as follows. In Sec. II, we will review the calculation of the azimuthal 
dependence of SIDIS up to twist-3 contributions in terms of the TMD parton distributions. In Sec. Ill, we extend 
the study of nuclear dependence of TMD parton distributions to the case of other TMD parton correlation functions, 
including twist-3. Using these results, we calculate the general nuclear dependence of (cos <fi). We will then illustrate 
the nuclear dependence with an ansatz of the TMD parton distributions in a Gaussian form and discuss the effect of 
transverse momentum smearing in the fragmentation. A summary will be given in Sec. IV. 

II. AZIMUTHAL ASYMMETRY IN SIDIS 

We consider semi-inclusive deep inelastic scattering e~ + A — >• e~ + q + X with unpolarized beam and nucleus or 
nucleon target. The differential cross section, 



o 



2 „2 



da = q V" v (l I') ^ - - - (7) 
can be expressed as a product of the leptonic tensor, 

i') = \\m' v + n /M - (/ • i')gH, (8) 
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and the hadronic tensor 
d 2 W^ 



d 2 k\ 



dk' 



■W^(q,p,k'); 



W^(q,p,k') = ±-J2(A\J„{0)\k',X)(k',X\U0)\A) 



x(27r) i S 4 (p + q-k' - Px ), 



(9) 



(10) 



where the superscript (si) denotes that it is for SIDIS. Here I and I' are the four momenta of the incoming and outgoing 
leptons, respectively, p is the four momentum per nuclcon of the incoming target A, q is the four momentum transfer, 
and k' is the four momentum of the outgoing quark. We neglect the masses and use the light-cone coordinates. The 
unit vectors are taken as, n = (1, 0, 0, 0), n = (0, 1, 0, 0), n±_\ = (0, 0, 1, 0), n±2 = (0, 0, 0, 1). We choose the coordinate 
system such that, p = p + n, q = —xsp + nQ 2 / (2xbP + ), and l± = |Zj.|nj_i, where xb = Q 2 /2p-q is the Bjorken variable 
and y = p ■ q/p ■ I. 

Since the azimuthal asymmetry (cos <p) in the kinematic region k± <C Q is a twist-3 effect, we need to calculate 
d 2 W flL , /d 2 k± up to the twist-3 level. The calculations have been carried out in Ref. [HI and we review it here for our 
later study of nuclear dependence. 

In the unpolarized SIDIS e~ + A — > e 
is given by, 



Twist -2 



d 2 k ± 

where the hard part is 
1 



q + X , the twist-2 contribution is independent of the direction of k± and 
H${x B )tf(x B ,k x ), (11) 



#(0)( 



XB 



Ap-q 



(12) 



with d^ v = — n^n" — n" n 11 as a projection tensor, and fg(x, k±) is TMD quark distribution function in the nucleus 
or nuclcon, 



dy d 2 y ± 
(2tt) 3 ' 



(A\i,(0)^-C(0,y)^(y)\A), 



(13) 



For brevity, we now work in the covariant gauge and use y to denote four vector of coordinates (0, y ,y±). The gauge 
link C(0,y) is given by, 

£(0,2/) = £H-w,d ± ;0,d ± )£(-^,y ± ;y-,y±); (14) 



-oo, y±;y 



,y ± ) = Pexplig f d£~A + {£- 



(15) 



Twist-3 contributions to the semi-inclusive hadronic tensor include 3 terms which can result in cos (f> azimuthal 
dependence, 



'd 2 W^~ 






+ 

Twist-3 


\d 2 wfr L) ~ 


+ 


\d 2 wfr R) ] 


. d2k ±- . 


Twist-3 


d 2 k ± 


d 2 k ± 


Twist-3 


d 2 k ± 



(16) 



Twist-3 



where the number in the superscript of W^ v denotes the number of gluon(s) involved in the multiple gluon scattering 
with respect to which the collinear expansion is carried out; the tilde above W denotes the results after collinear 
expansion; and the superscript L or R denotes the left or right cut respectively with respect to the initial gluon line. 
The contribution is given by, 



d 2 wfj 

d 2 k ± 



H$(x B )&V A (x B ,k ± )d>>°, 



(17) 



Twist-3 



where the hard part is given by, 

and the TMD parton correlation is defined as, 

^^)-J (27r) 3 £ 

After carrying out the trace, we obtain, 



.la 



(A\^(0)^-C(0,y)i>(y)\A). 



d 2 k 



x 



Twist -3 

where the twist-3 parton distribution /^(ib, k±) is defined as 



The contribution W^' R) and W^ L>I are the same in unpolarized SIDIS and are given by, 



d 2 w^ d 2 wfr L) 



d 2 k ± d 2 k ± Ap-q 1 P *> 

where w p p = g p — n p n p , Wjtv = lni>n p f^iv and the matrix element 

A {x,kx) = 



p+dy d 2 y ± eixp+y -. iS± .k ± (A |^( 0)£(0) y )D p (y)ijj{y)\A) , 



(27t) 3 

has two independent terms contributing to the hadronic tensor 



d 2 fc 



J Twist-3 



where ej_ P7 = e a p pl n a n^ . After carrying out the traces and make the Lorentz contraction with fcj 
respectively, we obtain, 



d 2 W$ 
d 2 k! 



J Twist-3 

where the parton correlation functions are given by 
kl^ )A {xM) = n a k p ± ^ A {x,ks_), 



p ■ q 



ikl^l )A {x, k ± ) = n a effk ±p <pV> A (x, k ± ), 



<p$ A (x,k x ) = 



^^e^ + ^-*-^(A|^(0)^£(0, 2/ ) J D p (y)^( 2 /)|A) 



<p$ A (x,k ± ) = J ^^e^ + ^-^-^(A|^(0)^£(0, 2/ )^ p (y)^(2/)|A) 
Equation of motion relates 

xf A ±(x,k ± ) = -(p ± L)A (x,kx) + tp ( l' A (x,kx), 
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so we have, 



~d 2 W^~ 




\d 2 w£) 


d 2 k ± 


Twist -3 


d 2 k ± 



d 2 kj_ 



Twist -3 



1 

P ■ 1 



[(g-u + 2x B p IJ )k^ u + (q v + 2x B p„)kj_ IJ ,]f£ L (x B ,k±), (31) 



Summing both twist- 2 and twist-3 contributions and contracting with lepton tensor L M „ , we obtain the differential 
cross section as, 



da 



-{[! + (!- y) 2 ]ff(x Bl k ± ) - 4(2 - y)^~V 



-Q- X BJq_L {XB,k±)C08<pj, 



dx B dyd 2 k± Q 2 y 

where the azimuthal angle <f> is defined by cos <f> = k± ■ l± . The azimuthal asymmetry at fixed k± is given by, 
2(2 - y)y/T=y \k±\ x B f^(x B , k±) 



(32) 



;cos <p) eA 



fMx B ,kx) 



i + (i-y) 2 k 

We can also calculate the transverse-momentum integrated asymmetry and obtain, 
2(2 - y)y/T= yj \ki_\d 2 k ± x B f^(x B , k±) 



((cOS(j))) eA = 



i + (i-y) 2 



f*(x B ) 



(33) 



(34) 



then, C = 1 , D a = d a , 
k±\/Q, which is just the 



where fq(x) = J d 2 k±f^(x,k±) is the usual quark distribution in a nucleon or nucleus. 

If we consider only "free partons " with intrinsic transverse momentum, i.e., setting g 
and xffj_(x,k x ) = ff(x,kx). In this case, (cos4>) eA = -2(2 - y)y/l - y/[l + (1 - y) 2 
result obtained in Ref. [33^ . 

We note that the above calculations apply to SIDIS of both nuclear and nucleon targets. All results are in the 
same form in terms of the TMD parton distributions inside a nucleus or a nucleon. The nuclear dependence of the 
azimuthal angle asymmetry will come from the nuclear dependence of the TMD parton distributions. In the SIDIS 
off a nuclear target, multiple gluon scattering between the struck quark and the target as contained in the gauge 
links can happen to different nucleons inside the nucleus. This will give rise to the nuclear dependence which we will 
discuss in the remainder of this paper. 



III. A-DEPENDENCE OF THE AZIMUTHAL ASYMMETRY 



Following the same approach in the discussion of the nuclear dependence of the twist- 2 TMD quark distribution 



fg(x, k±) in Ref. [241 ]. we can also express the general parton distribution, 



P + dy d 2 yj_ 



2- 

p+dy- 



(2tt)' 



- ik ^(A\i>(Q)^-C(Q,yW(y)\A), 



Jxp+y ^| ^_2L p Wx{y )-v kA 



i;(y-)\A)S^(k A 



(35) 



(36) 



2tt N 1 rw 2 

in terms of the collinear parton matrix elements involving the transport operator W±(y~) [Eq. ([2"]l]. where r a is any 
gamma matrix. Expanding the exponential term of the above matrix element in powers of the transport operator 
W±(y~) and assuming "maximum two-gluon correlation approximation" as in Ref. p4| . one can express the nuclear 
TMD parton distributions in terms of a Gaussian convolution of the same TMD distributions in a nucleon, 

1 (37) 



K{x,ki_) 



Aexp 
A 



™*Z(x,l ± ). (38) 

7tA 2 f J 

Note that in the derivation of the above result, one has considered the fact that matrix elements with odd powers of 
W±(y~) vanish for T a = 7+ while the matrix elements with even powers of W±(y~) vanish for T a = 7_i_. 

The nuclear broadening of the twist- 2 TMD quark distribution is a special case of the above general result with 
T a = 7 + . With the definition of the twist-3 parton distribution function f^_(x,k±) in Eqs. (|19j) and (|2~Tj). one can 
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multiply both sides of Eq. (|38j) by d ap and obtain, 

/£(*,*x) « / *ti£^e-$^r'*»f» ± {z,i ± ) (39) 

7rA 2F 7 fc^ 

'l + %jfc x • 9 fc .l / (40) 



Given the twist- 2 TMD quark distribution function f^{x,k±) and the twist-3 quark distribution f^ ± (x,£±) in 
a nucleon and using the above convolution, Eqs.(j31) and (|40p . one can then calculate the nuclear dependence of 
the azimuthal asymmetry (cos(f>) and ((cos<p)). In general, if the convoluted twist-3 TMD quark distribution is a 
decreasing function of the transverse momentum in the region of interest, the second term involving a derivative in 
the above equation will be negative. Therefore, one would expect the azimuthal asymmetry to decrease because of 
the multiple scattering in nuclei. 

To illustrate the nuclear dependence of the azimuthal asymmetry qualitatively, we consider an ansatz of the Gaussian 
distributions in k± for both the twist-2 and 3 TMD quark distributions, 

f?(x,k x ) = ±f?(x)e-*-' a , (41) 

f q N A^k ± ) = ±f?Ax)e- kl/ e. (42) 
The corresponding TMD distributions in nuclei are, 

fq(x,k±) » - A- X (x)e- k ^ a+ *>r\ (43) 
H ir(a + A 2 F) 



7r(/3 + A 2F y 



One can then calculate the the azimuthal asymmetry for SIDIS off both nucleon 

pJV 

l + (l-yf (3 Q f»(x B ) a(3 

and nuclear targets 



2(2 - y)y/T=jj a \k ± \ x B f^(x B ) r a-$ u2 \ .... 



/^^\ 2(2 - y)^^^ (3(a + A 2F ) |fc ± | *Wff>) r a - g 2 i 

The nuclear modification factor for the azimuthal asymmetry is then, 
(cos0) eA _/3 2 (a + A 2F ) r (a - p)A 2F (a + (3 + A 2F ) -> 2 



(cos<f>) eN a(/3 + A 2F ) 2 ' CXP l a/3(a + A 2F )(/3 + A 2F ) 
In the special case when a = ft, we have 



•I}. (47) 



{C °^ U = T+Jl^y) 2 Q f^xsT' (48) 



2(2-y)y/T=jj g |fc±| ^/glW 
l + (l-y) 2 /3 + A 2F Q 



(cos 4>) eA = T - Tt N , (49) 



(cos^ _ p m 



(cosip) eN j3 + A 



2F 
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Therefore, the azimuthal asymmetry (cos (j>) e A in deep inelastic eA scattering is suppressed compared to that in eiV 
scattering and the suppression is inversely proportional to the total transverse momentum broadening A 2 f- The 
suppression is independent of the transverse momentum k±. However, in general, the twist- 2 and twist-3 TMD quark 
distributions are not necessarily the same and their Gaussian ansatz might have different widths ft ^ a. The nuclear 
modification factor for the azimuthal asymmetry will then have non-trivial k± dependence. Shown in Fig. 1 are 
the nuclear modification factors for the azimuthal asymmetry when ft/a = 2 and 0.5, respectively, as functions of 
&2f/ol, at different transverse momentum k±. In the case ft > a, we see that the azimuthal asymmetry is suppressed 
and the suppression increases with the transverse momentum kj_. However, when ft < a, the suppression actually 
decreases with increasing k± and the azimuthal asymmetry could be enhanced for large enough transverse momentum 
kj_ . Therefore, the nuclear modification of the azimuthal asymmetry and its transverse momentum dependence is a 
very sensitive probe of the twist- 2 and twist-3 TMD quark distribution functions. 





FIG. 1: (color online) Ratio (cos <jf>) e A/(cos 0) e jv as a function of A2F for different k± and the relative width ft /a of twist-3 and 
2 TMD quark distributions. 

If we integrate over the magnitude of the transverse momentum k± , the averaged azimuthal asymmetry will only 
depend on the shape of the twist-3 TMD quark distributions, 

// j\\ 2{2~y)jT^ ^ft XBf^{x) 

((cos0>u = - ~2Q~fiw (51) 
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fN 

«cos0})eA = - T, m 7^J/ . (52) 



2(2-y)yT~i; Vtt/3 XBffAx) 
l + (l-y) 2 2QV£ + A 2F 



((cos</>)) eyt _ / P ,-„, 
((cos<«U ^ + A 2F ' 

We see again that ((cos </>)) eJ 4 is suppressed compared to ((cos (f>)) e N, and the suppression factor is inversely proportional 
to the square-root of the total transverse momentum broadening A 2 f- 

To take into account of the transverse momentum during the quark fragmentation process and its effect on the 
azimuthal asymmetry, we take another Gaussian smearing 

da eN^ehx = J da eN ^ eqX Dp h (z, k F ±_)dzd 2 k F (p h - zk' - k F1 _)d 3 p h , (54) 
for the fcfj_-dependence in the fragmentation function D q f* h (z, kp±), i.e., 

D F -> h {z,k FA _) = Dp h {z)— e - p *J a *. (55) 

ira F 

Consider the case that a = /3 and one has the azimuthal asymmetry for the hadron production cross section 



2(2-y)VT— y pz \p h± \ XBf?A*) 

{cosMeN = - i + d-yr ^ + a F —-j7W' (56) 



fnn A \ 2(2-y)VT— y pz \p h± | XBjfAx) 

{coscf) h } eA = — — (o . a — w~ n fNI \ ■ ( 57 

1 +(!-#) [P + A 2F )z z + a F Q f£(x) 

We compare the above results with those obtained without fragmentation, Eqs. (|48p and (|4^|. and see that we have a 
clear smearing effect on the azimuthal asymmetry in both e^iV and e~ A-scatterings. The smearing factors are given 
by 



(COS 4>h)eN 



(COS <p) eN 



(COS (ph)eA 



(COS <p) e A 



\p h± \= z \k ± \ /3z 2 + a F 

(/3 + A 2F )z 



\p h± \=z\k ± \ (fi + A 2F )z 2 + a F 
The suppression factor in eA compared to eN is given by, 

(COS 4>h)eA _ /3z 2 + a F 

(coscj) h } eN ((3 + A 2F )z 2 + a F ' 

After integrated over the magnitude of the transverse momentum, we have, 



(59) 



(60) 



((cos<fa)) eA = / Pz 2 + a F 
{{cOS(j)h}}eN \] W + A 2F )z 2 + a F ' 

IV. SUMMARY AND DISCUSSIONS 

Within the generalized factorization, we have calculated the SIDIS cross sections in terms of the TMD quark dis- 
tributions in a nucleon or nucleus up to twist-3. The azimuthal asymmetry (cos</>) in the small transverse momentum 
region depends on both twist-2 and 3 TMD quark distributions. By considering nuclear broadening of both twist-2 
and 3 TMD quark distributions due to multiple scattering between the struck quark and nucleons inside the nucleus, 
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we investigated the nuclear dependence of the azimuthal asymmetry. We found that the azimuthal asymmetry is 
suppressed by multiple parton scattering for most cases of the TMD quark distributions. The suppression is inversely 
proportional to the average squared transverse momentum broadening. The transverse momentum dependence of the 
suppression depends on the relative shape of the twist- 2 and 3 TDM quark distributions. Using a Gaussian ansatz, we 
found that the suppression factor decreases with the transverse momentum if the width of the twist-2 TMD distribu- 
tion is smaller than that of the twist-3 TMD distribution, while the suppression factor increases with the transverse 
momentum if the width of the twist-3 TMD distribution is smaller than that of the twist-2 TMD distribution. The 
suppression is independent of the transverse momentum if the twist-2 and 3 TMD distributions have the same width. 
Therefore, study of the nuclear dependence of the azimuthal asymmetry can shed light on the relative shape of the 
twist-2 and 3 TMD quark distributions. 
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ogy of Shandong Province and the Director, Ofhce of Energy Research, Office of High Energy and Nuclear Physics, 
Division of Nuclear Physics, of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. 
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